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Abstract —Two-way is a dominant mode of communication in 
wireless systems. Departing from the tradition to optimize each 
transmission direction separately, recent work has demonstrated 
that, for time-division duplex (TDD) systems, optimizing the 
schedule of the two transmission directions depending on traffic 
load and interference condition leads to performance gains. In 
this letter, a general network of multiple interfering two-way links 
is studied under the assumption of a balanced load in the two 
directions for each link. Using the notion of interference spin, we 
introduce an algebraic framework for the optimization of two- 
way scheduling, along with an efficient optimization algorithm 
that is based on the pruning of a properly defined topology 
graph and dynamic programming. Numerical results demon¬ 
strate multi-fold rate gains with respect to baseline solutions, 
especially for worst-case (5%-ile) rates. 

Index Terms—Two-way communication, scheduling, dynamic 
TDD, dynamic programming. 


I. Introduction 

Two-way is a dominant mode of communication in wireless 
systems, such as in uplink (UL)/downlink (DL) cellular com¬ 
munication. Although the two-way channel is the first known 
multi-user channel treated in information theory fT), the design 
of two-way links has traditionally been tackled by separating 
the problem into independent one-way DL and UL problems, 
respectively, that deal with rate adaptation, scheduling, etc. 

Nevertheless, recent works El-El have demonstrated that 
the presence of interference motivates the joint consideration 
of UL and DL, since the interference caused by a link on other 
links is different depending on the direction in which the link is 
active, i.e., on which node acts as the transmitter. Specifically, 
the referenced works have studied dynamic Time Division 
Duplex (TDD) for small-cell wireless systems. With dynamic 
TDD, each slot, of possibly different size in the frequency-time 
plane, can be assigned to either the UL or DL direction, de¬ 
pending on traffic load and interference conditions. References 
El and 01 outline the challenge of managing UL-DL cross¬ 
interference in the presence of dynamic TDD. The works (4), 
0 instead propose centralized and decentralized algorithms 
that optimize the switching time between UL and DL. Finally, 
0, 0 put forth heuristic solutions for the optimization of 
UL-DL slot allocation for fixed switching times. 

In this letter we consider a generic collection of two-way 
links, hence including, e.g., cellular and/or device-to-device 
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Fig. 1. (a) Example scenario with two macro-cells and one small cell, with 

direct channels (full lines) and interfering channels (dashed lines), (b) Example 
spin configuration in two subsequent frames, (c) Topology graph Q. 


(D2D) links, where the traffic load is balanced between the two 
transmission directions: a slot assigned to one transmission 
direction is always followed by one allocated to the opposite 
direction. This assumption is per definition suitable for appli¬ 
cations using interactive communication. An example setup is 
given in Fig. |Tja), that consists of three interfering two-way 
links operating over two-slot frames, with one slot allocated to 
each direction as per Fig. |TJb). We refer to the binary variable 
that determines the link direction in the first slot (the opposite 
to the one in the second slot) as the interference spin of a link 
(see precise definition in the next section). This allows us to 
propose an algebraic framework for the optimization of two- 
way scheduling, e.g., UL-DL slot allocation. The optimization 
consists of selecting the interference spins of all M links in 
order to maximize a general sum-utility of the signal-to-noise- 
plus-interference ratios (SINRs) of the two-way links. The 
optimal solution has a complexity that scales as 2 M . Besides 
the algebraic framework, our second contribution is a subop- 
timal, but efficient, algorithm that leverages the representation 
of the M interfering two-way links as a graph of M nodes. 
The proposed algorithm is based on pruning this topology 
graph according to a specific criterion and applying dynamic 
programming. The algorithm finds a spin configuration that is 
seen via numerical results to exhibit substantial performance 
gains over the baseline scheme in which the spins are not 
optimized. 
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II. System Model and Problem Definition 

We consider a set V of M synchronous two-way links, 
where the l— th link Vi consists of two half-duplex nodes L[ 
and Ri. The labeling of the two nodes in a link Vi as left, L[, 
and right, Ri, is arbitrary, see Fig. Q] Each link Vi uses TDD, 
and two-way communication takes place in a frame consisting 
of two slots. The first slot is used in one direction, either 
Li —» Ri or Ri —>• Li, and the successive slot is used in the 
opposite direction. Related to this, we define interference spin 
or, for short, spin of a link in a given frame: the 1-th link 
Vi is said to have a left spin or 0—spin if L; transmits in the 
odd slot and Ri transmits in the even slot; otherwise the link 
has a right spin or 1—spin in that frame. In the example in 
Fig. 0Jb), the link 14 has a left spin in the first frame and a 
right spin in the second frame. The spin of the Z—th link is 
denoted by si £ {0,1}. 

All devices use the same spectrum in a TDD manner, 
leading to interference among the concurrent link transmis¬ 
sions. Each device is backlogged with data, such that there is 
transmission in each slot and there is a continuous stream of 
slots. All links are slot-/frame-synchronous. Considering two 
links, say V k : L k — Rk and Vi : Li — Ri, we distinguish 
the following channel gains, see Fig. [2] (a) Direct channels: 
L k ->• R k , R k ->• L k , Li Ri and R t L u whose 
signal-to-noise ratios (SNRs) are defined as SNR^ fl , SNRf" L , 
SNR^ fl and SNRf L , respectively; (b) Interfering channels: 
Lk t Rh Rk ~^ L[, L k —> Li and R k —> Ri , whose 
interference-to-noise ratios (INRs) are defined as INR^/*, 
INR^ l , IN 111;/' and INR^ H , respectively. The interference 
from l to k, not depicted in Fig. [2] is represented in a similar 
manner. Even with perfect channel reciprocity, we generally 
have SNR k R ^ SNR£ l and INR^f ^INR^ L , since the 
nodes at the two ends of a link may use different powers. 
Some of these quantities can be zero: for instance, links may 
not interfere with each other, e.g. due to obstacles and yield 
INRjf =INR kl L = 0. 

In order to capture the interference-related features of the 
network, we define the undirected topology graph Q = (V, £), 
where the set of vertices V represents the links Vi, as defined 
above, and an edge exists in the edge set £ for a pair of 
links 14 and Vi in £ if and only if at least one of the 
interfering powers INR^f, INR^ L , INR^f, INR^, INR^f, 
INR^ l , INR^. l and IMl/E is non-zero. Therefore, an edge 
kl £ £ indicates that two links 14 and Vi interfere at least in 
one direction. Note that we identify edges via the indices of 
the connected links. The topology graph for the example in 
Fig. |TJa) is depicted in Fig. [Qc). Without loss of generality, 
the graph is assumed to be connected 0, since, otherwise, one 
could consider the different connected components separately. 

Denoting by © the XOR operation, the interference between 
two links 14 and Vi is fully specified by the relative spin r k f. 

r k i = r ik = s k © si. (1) 

The SINRs for any link Vi in the two directions Li —>■ Ri 
and Ri —» Lj can be written as a function solely of the 
relative spins of the interfering links, i.e., of r k i with kl £ £. 
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Fig. 2. Two interfering links Vfc : L*. — R & and V\ : Li — R{. The INRs 
are given only for the interference from Vi to Vk for clarify of illustration. 


Specifically, the SINR for the direction Li -A If is given as 

SINRf fl (r;) =- m - m - 

l + E fc : We£ [( 1 -^)INRf, il + r w INRjg R ] 

( 2 ) 

and for the direction Ri —> Li we have 

SiNR R/ 4r,l =_ SNR rl _ 

' l + E fe:Wef [(l-^)INR^+r w INR-]’ 

(3) 

where r; is the vector of spins r k i, kl £ £, for the link Vi. By 
rewriting the denominator of Q as 1 + kl ^ £ [INR,[ j ( r + 
r H {INR % r - INR^f)], we see that, if INR^ = INR^ fl , 
the relative spin r k i does not affect SINRf^rz). Similarly, if 
INR^ l = INR^ l , then r k i does not affect SINRf L (ri). 

Given a topology graph Q , definition ([[} implies that the 
vector r = [iq • • • r M ] of all relative spins satisfies the 
following properties: (Cl) Symmetry: r k i = ri k \ (C2) Par¬ 
ity check on cycles: For any cycle in the topology graph 
III 2 , hh---, l N h G £, we have the parity check equation 

rip 2 © n 2 i 3 © ... © ri N _ 1 i N © n Nh = 0 . ( 4 ) 

It can be easily shown that Cl and C2 are also sufficient to 
guarantee the existence of a vector of spins that satisfies 0 
and hence they characterize the set of all relative spin vectors. 
A simple consequence of Cl and C2 is that the specification 
of the relative spins on the edges of any spanning tree 0 on 
the topology graph Q is sufficient to obtain the entire vector r. 
Every edge kl either belongs to the tree and hence r k i is fixed, 
or, if not, then edge kl forms a unique cycle with a subset of 
the edges of the tree and the relative spin follows from IQ}. 

We are interested in finding the spin configuration that 
maximizes a sum-utility function of the SINRs across all links. 
The problem of interest is hence formulated as 

M 

max V Ui(SWRjf R (ri), SINRf^iy)) s.t. Cl and C2 

r ^ 

1=1 

(5) 

where the maximization over the relative spins is subject to 
the constraints Cl and C2, and we have fixed non-negative 
utility functions for each link V). The utility function 

(/;(•,•) i s assumed to be non-decreasing in the arguments and 
will be written as Ui (ri ) in order to simplify the notation. An 
example is the two-way sum-rate, which is given as Ui(vi) = 
log(l + SINRf^r 1 )) + log(l + SINRf L (r z )). In principle, a 
link can change the spin per frame; however, in Sec. HY1 we 
will treat the case in which a link spin is determined based on 
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SINR variables that include long-term fading and update the 
the scheduling decisions on a large time scale. 

III. Optimizing the Interference Spins 

The maximization in (0, in principle, can be carried out by 
exploring all possible configurations of relative spins. Due to 
the constraint (|4j, the complexity of exhaustive search scales as 
2 M_1 , becoming impractical as M increases. Here we propose 
a suboptimal, but computationally more efficient optimization 
strategy that consists of two steps: 1) Construction of a 
maximum relative-interference spanning tree T over the graph 
Q\ 2) Dynamic programming spin optimization on T. 


with the average of the interference powers that would be 
observed if r^i = 0 or = 1. This is justified by the fact 
that, in light of the choice of the weights ([6]), it is expected 
that the two values INR^ and INR^/* are similar. The same 
approach is used to approximate @, leading to the approx¬ 
imation SINRf i (r Ci i,np i ). For ease of notation, we write 
Ui{r Cl i,ri Pl ) = Ui (SINRf(r Ci i,n Pl ), SINRf L (r Ci ;, )). 

We now aim at optimizing problem (0 with utilities Ui in 
lieu of Ui for all l £ V. The proposed dynamic programming 
solution starts from the leaf vertices and proceeds according 
to the (partial) order defined by the tree until the root is 
reached. In particular, each vertex V/ calculates the message 
p l = {p{,p\) for its parent vertex Pi, where 


1) Maximum Relative-Interference Spanning Tree 'T 
For each edge kl in Q , a weight Wki is assigned as 

wm = max(|INR^ - INR“|, |INR^ - INR“|, 

|INR« fl - INRf/|, |INRf fe L - INRf/'l). (6) 

This weight evaluates the maximum change in interference 
powers that is affected by the selection of the spin rri ■ A 
maximum spanning tree T is then constructed with respect 
to the weights in ((6)), retaining only the edges that have the 
largest contributions to the relative interference powers and 
pruning the remaining edges; e.g., for the example in Fig.[TJc), 
one of the three edges will be removed based on the criterion 
©, thus obtaining T with two edges. We refer to T as the 
maximum relative-interference spanning tree. The tree T can 
be constructed in a centralized or distributed way GO M, but 
the implementation details are out of the scope of this paper. 
Finally, a root vertex is arbitrarily selected so as to make T a 
rooted tree. Note that each vertex Vi in T has a single parent 
vertex Pi (except for the root vertex, which has no parent), 
i.e.. Pi is the only vertex in T along the unique path to the 
root. Vi can be the parent for multiple children vertices C/. 


for i = 0,1. The maximization in © is over the relative spins 
corresponding to edges stemming from the child vertices Ci. 
We denote a solution of the problem (0 as r' Cil for i = 0,1. 
Note that the leaves have no child vertices, and hence, for 
every leaf V), the message is calculated as p\ = Ui(ri Pl = i) 
for i = 0,1. Instead, the root vertex t, which has no 
parent, solves the problem max rCtt ((%(r Ctt ) + Y.keC t Mfc**) 
and obtains an optimal solution r c t t- A complete solution f is 
finally obtained by backpropagation: starting with the children 
of the root vertex, each child vertex k £ Ci of a vertex V) 
selects the solution until the process reaches the leaves. 

The complexity of the maximum spanning tree construction 
scales as GlogM J8] J9), where G < M(M — l)/2 and 
M are the number of edges and vertices in Q, respectively. 
The order of complexity of the relative-spin optimization is 
M ■ 2 d , where D is the maximum number of children of a 
vertex in T. While in the worst case, D = M — 1, for a 
typical topology graph Q , we found D to be much smaller 
than M, leading to significant complexity saving with respect 
to exhaustive search, see Section m 


2) Dynamic Programming 

Having constructed the spanning tree T, we now proceed 
to optimize only the relative spins corresponding to the edges 
kl £ T via dynamic programming. As discussed in Sec. [Til 
the spins for all the remaining edges can then be immediately 
calculated via (|4j. In order to allow the optimization to be 
limited only to the edges in T, we approximate the SINRs in 
(O and 0 so that they only depend on the relative spins for 
the edges in T- To this end, denote by r cp the vector of the 
relative spins for the edges connecting the child nodes C/ with 
Vi and, similarly, by ri Pl the relative spin between Vi and its 
parent Pi. The SINR (0 for link Vi is approximated as 

SINRf fl (r Cl i,ri Pl ) = - - SNR; - , where (7) 

1 + Z^fc: kl e £ IbJRfcJ 

(1 - r fe z)INR^ + r w INR« fl k £ C t U {Pi}, 
INR^ + (INR^f + INR^ fl )/2 else. 

( 8 ) 

In other words, the interference contribution for the edges 
kl £ £ \ T that do not belong to the tree is approximated 


IV. Numerical Results and Concluding Remarks 

For performance evaluation, a small-cell set-up is consid¬ 
ered over a 100m x 100m area. M links, or equivalently 2 M 
nodes, are generated as follows: M transceivers are placed 
uniformly in this area and chosen equiprobably to be either 
an L— or R— node; then, the opposite node of the link (R 
or L, respectively) is placed at a uniformly selected random 
angle with distance d. Note that the opposite node may lie 
outside the area at hand. Two types of links are considered: 
(i) symmetric D2D links, with d = d s = 10m; (ii) asymmetric 
femtocell links with d = d a = 50m. These links differ as 
explained below. The long-term SNRs and INRs that are used 
by the algorithm to determine the configuration of relative 
spins account only for large-scale fading and path loss. The 
long-term SNR parameters for a symmetric link V/ are given 
as SNH / /i ’ l = SNRf= SNH. s /l;. The transmit power is only 
adapted to the path loss in order to get SNR s = 20dB, while 
pi represent independent log-normal shadowing with standard 
deviation 8dB. For an asymmetric link V'/, the SNRs are given 
as SNRf fl = SNR**# and SNRf L = SNRf^ with 
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Fig. 3. 5%-ile two-way sum-rate (for a bandwidth of 10MHz) versus the 
average two-way sum-rate for different number of links M. 


SNR“ = 20dB and SNRf L = lOdB due to the fact that the 
power of the femto-base station (L) is different from the power 
of the mobile user ( R ). The long-term interference power 
caused by the node Xi to the node Y/ c , where X, Y £ {L, II}, 
is given as INR^ r = SNR ^ [d z /) v , where z = s 
if Vi is symmetric and z = a otherwise, XX £ {LR , RL}, 
r/ = 4 is the path loss exponent, dj} : Y is the distance between 
node Xi and Yk, and accounts for log-normal shadowing 
with standard deviation 8dB. Note that = /3j^ A . 

For a fixed spatial configuration and large-scale fading, the 
spins are optimized with the proportional fairness utility across 
the links, so that Ui(ri) = log(log(l + SINRf fl (r;)) + log(l + 
SINR^ L (rj))). The proposed algorithm is referred to as MST- 
DP (Maximum Spanning Tree-Dynamic Programming). One 
reference, when computationally feasible, is the the optimal 
performance based on exhaustive search; the other is the 
average performance with uniform random spins. For a fixed 
spatial configuration, large-scale fading and relative spins, we 
evaluate for each link Vi the instantaneous two-way sum- 
rate log(l + sinrf^r;)) + log(l + sinrf L (r;)), where the 
instantaneous SINRs sinrf fl (rj) and sinrf L (r;) are calculated 
via (0 and 0 by multiplying the corresponding long-term 
SNR and INR variables by unit-power Rayleigh variables, 
independently drawn for each par of nodes in each frame. 

In Fig. [3] we plot the 5%-ile of the two-way sum-rate, 
which is a common measure of worst-case rate, against the 
average two-way sum-rate, for different number M of sym¬ 
metric links and a bandwidth of 10MHz. Note that having 
more links M enables the spatial reuse, and hence the sum- 
rate, to be improved, but, at the same time, it creates more 
interference that reduces the worst-case rates. It is seen that 
spin optimization can significantly outperform random-spin 
assignment, especially in terms of 5%-ile rate, which, as an 
example, is nearly doubled for M = 10. 

Fig. a elaborates further on the 5%-ile performance gains 
attained by optimized spin assignment with respect to random- 
spin selection as a function of M. Note that the high number 
of links is justified by the ultra-dense wireless scenarios 
envisioned in the upcoming 5G wireless systems Q. Two 
curves are shown, one for all symmetric and one for all 



Fig. 4. Gain in 5%-ile two-way sum-rate for the proposed MST-DP algorithm 
with respect to a random spin selection versus the number of links. 

asymmetric links, respectively. For computational feasibility, 
only the MST-DP algorithm is considered. It is seen that, 
while the performance gains are substantial in both scenarios, 
optimizing the interference spins is particularly advantageous 
for asymmetric links. This is expected, since controlling the 
interference caused by the nodes that transmit with larger 
power has a more pronounced impact on the worst-case two- 
way sum-rate. Finally, the MST-DP algorithm, while inferior 
to the exhaustive search, still provides significant gains over 
random spins. We observe that, in terms of complexity, even 
with M = 100 nodes, the value of D was found to be always 
D < 9 and to be, on average less than 6. 

The main conclusion is that the degree of freedom offered 
by the interference spins leads to remarkable throughput gains, 
especially for the worst-case (5%-ile) rates. 
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